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Abstract 

We introduce graph state basis diagonalization to calcu- 
late the coherent information of a quantum code passing 
through a Pauli channel. The scheme is 5000 times faster 
than the best known one for some concatenated repe- 
tition codes, providing us a practical constructive way 
of approaching the quantum capacity of a Pauli chan- 
nel. The calculation of the coherent information of non- 
additive quantum code can also be greatly simplified in 
graph state basis. 

PACS number(s): 03.67.-a, 03.65.Ud, 
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The quantum coding theorem for noisy channels [T] 
[2 [3] [3] states that the quantum capacity Q{Af) of a 
channel Af is given by regularized coherent information: 
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Q{N)= lim -max/(p„,A/'^'' 



(1) 



the r.h.s term has long been known an upper bound for 
Q{N)^ which is the content of the converse coding the- 
orem [2] [3] . The direct coding theorem, stating that 
r.h.s tem is actually attainable, has been strictly proven 
by Devetak Q. The coherent information I{p,J\f) of a 
state p with respect to the noise Af is defined by 



I{p,Af) = S{Af{p)) - S{lA<E>Af{'i'Ap)), 



(2) 



where S{g) = —trg\og2g is the von Neumann entropy, 
^Ap is a purification of p, and Ia is the identity operation 
on the ancilla system A. The last term, S'(Ja'8'A/'(^Ap)), 
is the entropy exchange Se {p, Af) of p with respect to Af. 
For a code state p'^ of n qubits, denote a'^ = 
A/'«"(p^),a^^ = Ia (E) Af'^"-i\'^^^){'i'^^\), where 
1^'^'-^) is the purification of p*^, the coherent informa- 
tion of the state p*^ with respect to the noise A/"**" per 
qubit will be I^^ = i(5(o-C^) - 5(a^^)). Thus /^^ is 
the lower bound of Q{Af) according to quantum noisy 
coding theorem. It is known that the one-shot capac- 
ity Qi(Af) = maxp I{p,Af) is exactly the maximum rate 
achievable with a non-degenerate code for Pauli chan- 
nel Af. That Q{Af) > Qi{Af) is then established by the 
construction of a massively degenerate code, this was ac- 
complished in the work of [5] [B] for depolarizing channel 
and [7] [5] for some Pauli channels. It is not known which 
quantum code achieves the quantum capacity for a chan- 
nel that is neither degradable nor anti-degradable. Pauli 
channel with proper channel parameters is an example of 
such channels. So we need to check all possible codings to 



seek the maximal coherent information, this is an awful 
work in the viewpoint of just working out the quantum 
capacity. However, the history of classical communica- 
tion tells us that coding is the really important thing 
even when the capacity is known. The aim of this paper 
is to provide a scheme to work out the coherent informa- 
tion for a quantum code with respect to a Pauli channel. 
Graph state basis. — A graph G — {V^ F) is composed 
of a set V oin vertices and a set of edges specified by the 
adjacency matrix F, which is an n x n symmetric matrix 
with vanishing diagonal entries and Fat, = 1 if vertices 
a, h are connected and V ah — otherwise. The neighbor- 
hood of a vertex a is denoted \>y Na= {v &V |Fa« = 1 }, 
i.e, the set of all the vertices that are connected to a. 
Graph states [5] [IH] are useful multipartite entangled 
states that are essential resources for the one-way com- 
puting [11] and can be experimentally demonstrated [12] . 
To associate the graph state to the underlying graph, 
we assign each vertex with a qubit, each edge represents 
the interaction between the corresponding two qubits. 
More physically, the interaction may be Ising interaction 
of spin qubits. Let us denote the Pauli matrices at the 
qubit a by Xa ,Ya, Za and identity by la . The graph state 
related to graph G is defined as 



\G) = n u-b 
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^(-l)^^^'^"!^) (3) 
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where |/i) is the joint eigenstate of Pauli operators Za 
(a € V) with eigenvalues (—1)^", |+) is the joint +1 
eigenstate of Pauli operators Xa { a £ V) , and Uab 
{Uab = diag{l, 1, 1, —1} in the Z basis) is the controlled 
phase gate between qubits a and b. Graph state can also 
be viewed as the result of successively performing 2-qubit 
Control-Z operations Uat to the initially unconnected n 
qubit state l-l-)^, . It can be shown that graph state is the 
joint -|-1 eigenstate of the n vertices stabilizers 



Ka=Xa Yl Zb-.^XaZN^, a £ V. 



(4) 
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Meanwhile, the graph state basis are \Gki,k2,---kr,) 



Ua^y Z^^ \G) , with ka 
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ki,k2 
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0, 1. Thus 

= (-1)'° \Gk,M, 



(5) 



Output state in graph state basis. — Consider the input 
code state p'-^which is diagonal in graph state basis, that 
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7Tk\Gk){Gk\ 
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with TTfc the probabihty of state the graph state \Gk) ,0 < 
TTk < 1 and '^i^T^k = 1, we have denoted \Gki.k2,---k„) — 
\Gk) with the conventions oi k = J^i ki'^"~^ and 
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J2{h(sj.)r 



(7) 



The purification state can be Ivl/'^'^) = 

Sfc V^l^fc)^ l^fc) ■ ^^ Krauss representation, Pauh 
channel map Af acting on qubit state p can be written 
as J\f{p) = fp + VxXpX + PyYpY + PzZpZ, where 
Px(y,z) G [0, 1] are the probabihties, f — I — Px — Py — Pz 
€ [0,1] is the fidehty of the channel, and X,Y,Z 
are the Pauli operators. For depolarizing channel, 
Px — Py = Pz = P, f = ^ — 3p. For n use of channels 
with n qubits input state p„, we have the output state 

A^^"(Pn) = EaVaEaPnEl with TJa = p-'-'-' T^^PijP'z 

for Ea — X^Y^ZK Then the joint output state of 
p^ and the ancilla is ct"*'^ = J2ij V^^i^il^i) Ai^jlA^ 
®AA«"(|G.)(G,|) = E.,V^f^Ea^a(|G.)^(G,|^) 
(^{Ea\Gi) {Gj\El). In graph state basis, we 



have 



^AC 



(Gs|^ (G„ 



k^^lG*)^ 



G, 



E., V^^T.a Va {Gm\ E, |G.)_(G,| i?t \Gi) 5 ^Ay ThuS 

the matrix elements of the joint output state in graph 
state basis are 

^t^S = ,/^u^,Y.'^a{G.n\Ea\G.){G,\El\Gi) 

a 

= ./^UT^Y.'^a{G\Z^'^'iEaZ^^\G) 

a 

X (G| z(^)i?lzW \G) , (8) 



where we have denoted Z'^^'' — JJcev ^c" 
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cording to the orthogonality of graph state basis, 
(G| Z^-^^EaZ^^ \G) = except Z^-^^EaZ^"^ = K, up to a 
factor of ±1, ±i (the factor will be omit hereafter because 
it has no effect in the elements of cr"^'^) , for some Ks £ K 
(the group with generators of all Kc, c G V, an Abelian 
group, the vertices stabilizer group), . Thus we have 
Ea = Z^"^'^KsZ^\ so Z(j)£;tzW = Z(j)zWXsZ(™)zW. 
It can be written as Z^''®^'>KsZ^"^®^^ by our convention 
([7]). Since Z^^®^^ either commutates or anti-commutatcs 
with Ks, we have Z^^'>ElZ^^'> = iX^Z^*®^'®'"®'). For a 
non-zero (G| Z^^'>ElZ^''> \G) , we should get 



i®j ®m®l^O. 



(9) 



Let k = i(B m, then l = i(Bj®m — jQ) fc,and m — i(B k 
so the possible non-zero elements are in the form of 



AC 

^tii(Bk),JiJ(Bk) 
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a|B„ez<'=)A: 



where Pa = with E^ = Z^'+'^'iC^Zf') such that 
2'('®J) communicates with i^s,and Pa — \ with Ea = 
Z'-*+''^KsZ'-''huch that Z"^^®^^ anti-communicates with 
iiTs. Notice that the joint output state can be block diag- 
onalized according to fc. So that what left is to diagonal- 
ing each block with give k. The routine of calculating the 
non-zero elements is like this: (1) to list all the elements 



of K, (2) to determine P according to the commutator 
of Z(*®-')and i^,(3) to muhiply Z^''^ to obtain the coset 
E = Z^'''>K of the Pauli group with respect to its sub- 
group K and to determine rj, 

Notice that for i' ® i' — i ® ?', we have 
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i(i(Sk),](](Bk) 



/^Tf-TfJ 
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^AC 



t'(i'(Sk),]'{]'<Sk) 



/.jThW]- This 
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property is very useful in further diagonalizing the sub- 
matrix for stabilizer code of equal probability. 

Stabilizer code. — In a graphical quantum error- 
correction code, each codeword can be written as \Gk) = 
ZW|G) = ricev^MG)- To encode is to properly 
choose some of the \Gk) in order to form the code. A 
code is thus completely characterized by the set of k for 
a underneath given graph. For stabilizer code encod- 
ing q qubits into n qubits, all the 2'' chosen Z'^^^ forms 
a. group, each Z^'''^ is self inverse. Without loss of gener- 
ality, we use the binary vector (fci, ^2, . . . , fcg) to charac- 
terize the codeword of the stabilizer code. Then quan- 
tum stabilizer encoding is an encoding of classical binary 

' j if 
1 j' = i (B j, when i,j, i' ,j' correspond to codewords. 



serial (fci, fc2, . . . , kq) into binary serial (fci, ^2, . 
Denote k = Ez'=i^i2'~^, we have i' (B j' = 



A matrix M with M- = Mj, 



)i,j®i' 



can be diagonal- 



ized with Hadamard matrix Hq. Its eigenvalues will 
be yf2fl{HqM)jQ. Each block of a^^ can be written in 
the form of M for a code with a priori equal probability 
■Kj = 2""?. Hence, the eigenvalues of each bloc k of a^'^ 
will be 



At 



1 



{HqM)j-,. 



(11) 



The channel output state is simply a^ = TrA<J^'^, the 
matrix element of a in graph state basis is 



km 



EAC c 

'^i{i(Sk),i(i®kYkr: 



(12) 



In graph state basis, the output state of a Pauli channel 
with diagonal input is still diagonal, as mentioned in |13) . 
Thus far, we have obtained all the eigenvalues for calcu- 
lating the coherent information of the input of a priori 
uniform distributed stabilizer code with respect to Pauli 
channel. 

Concatenated repetition codes in depolarizing 
channel. — Depolarizing channel is a special case of 
Pauli channel, we will mainly deal with the depolarizing 
channel, the results can easily be extended to generic 
Pauli channel. One way to show Q{N) > Qi{M) is 
to find Q{Af) > for very noisy channel Af where 
Qi{N) = 0. Some codes that were shown to allow 
correction in the range of Qi{N) = consist of an 
ni qubit bit flip code concatenated with an 712 qubit 
phase flip code [7] [B] [5]. These have been called "ni 
in 712" codes [5], since each of the 712 blocks of the phase 
flip code consists of an ni qubit bit flip code. One of the 
examples is the famous Shor [[9,1,3]] code which is the 
"3 in 3" code. The codewords of "ni in 712" code can be 



obtained similarly as in Shor [[9,1,3]] code, they are 
1 
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.('|q'8i"i\ I h '8"^l\^(5 



(13) 



where 



L=(|0®"i) ± |l®"i)) are the ni-partite GHZ 
states. The "rii in ri2" codeword is the repetition of GHZ 
state. Comparing with the definition ([3]) of graph state, 
we find that the underneath graph for the code can be 
the "forest" graph. The graph contains 712 independent 
and identical subgraphs, each subgraph of n\ vertex has a 
tree structure with the root vertex connecting with all the 
other vertices and no other links exist. The vertices will 
be numbered in the following order: the j\ ~ th leaf of the 
J2—th tree is numbered as the {J2ini — l)+ji) — th vertex, 
the root of J2 — th tree is numbered as the j2'T'i — th ver- 
tex. We have |0) = \G) and |T) = Z„,Z2„i • • • Z„,„, \G) 
, where \G) is the graph state of the "forest" graph. The 
input state to the channel is p*" = ^(|0) (0| + |l) (l|), 
where the equal a priori probabilities are assumed. 

To simplify the analysis, we consider the case of 712 = 1, 
the tree graph first. Now |0) = \G) and |T) = Zn^ \G) . 
The vertices stabilizer group K can be divided into 
its subgroup if with generators Ki^ . . . ,Knj^-i and 
the coset K° — Kn^K'^, so that all the elements of 
K'^ commutate with Z^ and all the elements of K° 
anti-commutate with Z„j . The index i and j can be 
(0, • • • ,0, 0) or (0, • • • , 0, 1). Correspondingly, i and j can 
be or 1. In the basis of |0) and |l) ,each block of the 
joint output state a^'^ is a 2 x 2 matrix in the form of 



Ve + Vo Ve - Vo 
Ve -Vo r]e+ Vo 



(15) 



where r]e = X^a '^a with the condition of Ea G Z^^^'K'^; 
r]o = J2a''la with the condition of Ea £ Z'^'''>K° for the 
k~th block. For each block, rje and ijo are the eigenvalues 
ofa^^. 

For the fc = block of the joint output density matrix 
in Eq. (ITU)) , we have E = K. Now K*^ is generated by 
Ki = XiZni {I — 1, • • • ,rii — 1), the group elements of 
K'^ can be classified as two classes: the even class and the 
odd class. The even class have even generators, the group 
elements will be I, X1X2, • • • X„^_2X„j_i, • • • , the con- 
tribution to ?7e is /"I + ("'2~^)/"'"V+ ("V^)/"'"V + 
•■■ = ^[(/ + P)"'"^ + (/ - P)"'"^]/- The odd class 
have odd generators, the group elements will be 
XiZni,X2Zni, • • ■ , Xni-iZni, XiX2X3Zn^ • • • , the Con- 
tribution to 77e is ("VV"'"V +("VV"'"V + ■■■ 
= h[if+P)'''''-if-pT''^]P-^<^'^otex = p/f = p/il~ 
3p),wehaver]e{k^Q) = ^[(1 + j:)"i -h(l-a;)"i]. While 
for rjo, since Km — Xn^ Z1Z2 • • • ^„^-_i,the number of the 
Pauli error operators in each element of K° ~ K^K'^ is 
ni, so r]o{k = 0) = 



-ip"! 



niOTii — 1 rii 



r^2 



The next block we should consider is with the coset 
E = Ziii:,with E'' = ZiK" and E° = Z^K° . The 
group if now can be classified as two class: the 
class without K\ and the class with K\. The class 



without K\ contributes to r\e with ^^[(1 + a;)"^^"'^ + 
(1 — .T)"i~^]a;,the last factor x comes from Z\ in 
E'^ = Z\K^ . The class with K\ contributes to yye 
with -^^[(1 + x) 



.x)"i ^ — (1 — a;)"i ^\x, where the set 



ZxKx{I,X2X-i, • • • } contributes ^[(1 + x^^' 



(1 



xY^- 



and the set Z\K\{X2Zn^^X'iZn^^- ■ -^ con- 



tributes ^[(l-|-x)"i-2 
r;e(fc = 2"i-i) = /"i(l 



(1 — a;)"i '^\x. Thus we have 
X. While for ?7o, we have 



\ni-l 



E° = ZiK° = ZliC„,i^^ The class without Ki con- 
tributes to rjo with /"i2"i^^a;"i^^,the class with Ki con- 
tributes to 770 with /"i2"i-2a;"i, hence rio{k = 2"i-i) 
= /"i2"i-2(i + x)x'^^-^. For all the cosets E = ZiK 
(I — 2, . . . , rti — 1) , we obtain the same results of ?7e and r\o 
as in the case of the coset Z\K. A further research shows 
that for coset Zi^Zi^ ■ ■ ■ Zj^^iC with 1 < li < ni — 1, we can 



get rje = /"i2'"-i(l -|-x)"i-™x™,7]o = /"i2"i 
x) 



-Hi 



7n^ni —m 



For the coset E 
-^^[(1 +2;)"i - (1 



Znj^K, we can get rjeik — 1) 



1) 



pii2"i 



-l^ni 



^Zn-^K with 1 < li 



< 



r^lVoik 
Similarly, for coset Zi-^ Zi^ ■ ■ ■ Zi^ 
ni - 1, we have r/e(fc = 1) = tm+i,r]o{k = 1) = i„i-m+i. 
The eigenvalues of the joint output state a^'^ of tree 
graph case are 

to = ^[{i+xr + {i-xr], (16) 

h = ^[(l + a;)"^-(l-a:)"^], (17) 

t„,+i = ri2™-i(l + a;)"^-™x", l<m<ni(18) 

The eigenvalues of a^'^' whose value is t,„+i(l < m < 
ni - 1) has the degeneracy 2{"^~^) = 2C™_i. From Eq. 
p^ . the eigenvalues of the output state a'-'^are ^{rie{k) + 
Voik)) +i(77e(fc + 1) + Vo{k + 1)), which are i(io + ti + 
2i„i+i) and i(im+i -|- tm-m+i) for (1 < m < ni - 1). 

Denote the "forest" vertices stabilizer group as /C,with 
its subgroup ICj^ for J2 — th tree. Let's split the "forest" 
vertices stabilizer group /C into two parts according to the 
commutators of the elements and Zni^2ni ■ ■ ■ Z„2ni ■ The 
commutator of the element of /C and Znj^Z2ni ■ • • Zn^m 
for "forest" graph can be reduced to the product the 
commutators of the corresponding piece of JCj^oi and 
Zj2ni- When the element of JC commutates with 
ZniZ2ni • • • Zn^nx j the number of the trees with anti- 
commutator of its corresponding section of the element 
of Kj^ and Z^^n^ should be even. When the element of 
/C anti-commutates with Zn^Z2nx • • • ^n^nx^ the number 
of that should be odd. For a given coset E = Z^^^K. of 
the "forest", the "forest" coset head Z^^") is composed 
of ri2 sections, each section is the coset head of the tree 
graph case. Since the trees are identical, the types of 
the sections can be denoted with i, Z^ , Z\, Z^Z^ , ■Z^i^2, 

■^1-^2-^ni, • ■ ■ , ■^1-^2 • • ■ -^ni-l, ^1^2 ' ' ' ^ni-l-^ni ■ Here 

type Z\Z2- ■ ■ Zrn (1 < fn < ni — 1) represents all the 
cosets whose Z operator numbers are m. While the eigen- 
values of the cosets Z\Z2- • ■ Z^a and Z\Z2- ■ ■ Z^Zn^ are 
equal, we can further simplify the types of coset head 
sections as 7, Z^ , Zi, Z1Z2, • • • , Z1Z2 • • • Zm-i, with de- 
generacies 1,1, 2Cm-i, 2C^j_i • • • 2C"lZl, respectively. 
For a given coset head Z, suppose the number of trees 



with coset head / is /ojthc number of trees with coset 
head type Z^^ is ^i,the number of trees with coset head 
type Z\Z2 • ■ ■ Zm is Im+x- The total number of the trees 
is 7i2 — X]m=-i ^m+i- For a particular element of /C, con- 
sider the type Z1Z2 • • • Z^ trees, suppose there be Sm+i 
of the trees with their corresponding piece of /Cjj anti- 
commutating with their own Zj^m, this element of /C 
should contribute to the eigenvalues of the joint output 
state of the concatenated code with 



Table 1 



ni — 1 



Jo-SoJl-Si.So + Si TT .im + l-Sm+l, 

''0 '•1 '•rii + 1 11 '•m+1 '' 



+ lj.Sm+l 



ni— 771+1 



(19) 



Summing upon all elements of /C, we arrive at 
V" + V° - (to+t„,+i)'"(ti+i„i+i)'^ 



771 — 1 



_[_[ (^771+1 



t 



Til— m+1/ 



(20) 



7]^ - T]" 



Denote ao± 



m— 1 



tni + l) ^ 



771 — 1 

X Yl (trn+l 

777—1 
to ± ^771 + 1, 0-l± = 






(21) 



^1 ± ^771 + 1, a(7n+l)± 



777+1 



±t 



771-777 + 1 



(2 < ?7i < rii — 1), notice that a 



(777+l)± 



±a(„j_m+i)± for 2 < TO < rii — 1, we may only use Um+i 
with 1 < m < [ "^^"'" ] to specify the eigenvalues be- 
sides ao±, ai±. Consider the factors {tm+i+tm-m+iY"'*^ 
and {tn^-m+i + i7n+i)'"i""+' in Eq. (gOl), their de- 
generacies are (2C^^_i)'"+i and (2C"^'l™)'"i-'"+i, re- 
spectively. For 1 < m < [ "^^"'" J , let ^7„+i = 
Im+i + ^771-777+1, the degeneracy of a(„+i)+ should 

be Err;iUCM':;\(2c™-i)''"+^(2C-;")^'"+-'-+^ = 

(2C™_i -h2C;;^'l™)^'"+i = (2C™)^'"+i . For even m and 
TO = ^,let ^777+1 = ^777+1 the degeneracy of a(„_,_i)+ 
= 2t™+iwin be (2C™_i)'^-+i = (C™)^"+i. Thus the 
number of coset that gives the same 77*^ + 77° is 

r(»i-i)/2i 

d(M)=2'^n2! n T— T(^n:)^"^' 



M777+1 



(22) 



where vector fi = (^0, A^i, ■ • • , ^[{771-2)72]), h = 

Y-L(ni-1)/2J „ ■.^^ ^r(77i-l)/2] _ , 

2^777=1 A*777 + l, Wltll 2^„^_i A^777 + l " "-2 aUCl 

A'o.i = lo.i- 

Let 7i3 = [(ni — 2)/2] + 1, the eigenvalues of the joint 
output state a^'^ of "rii in 712" code can be written as 



-i "-3 't3 

777 — 777—0 

with degeneracy d(/z). According to Eq.([T 
values of the output state a'-' are 

773 



(23) 
the eigen- 



^'(m)-^«^ 



n 



'l + ^O+y 11 "-m+i 
777=2 



(24) 



also with degeneracy d (fi) . The coherent information of 
"7 in 712" code per channel use should be 
1 



tCA/ 



7ll7l2 



E 



d{^^)[-v' {fJ^)^os2v' {^^) 



'V+ (m) log2 V+ (/^) + V- (m) log2 V- (m)(25) 



Til 


?^2 


3Pmax 


9 


334 


0.19080163947 


10 


113 


0.19005842449 


11 


828 


0.19047624366 


12 


309 


0.18981106971 


14 


812 


0.18953710664 



Examples. — Consider "7 in 7^2" code, we have d{p) — 
_^ 14/^2 42^70/^4^ ao± = i[(l - 2p)7 + (1 - 4p)7 ± 

(2p)7],ai± = i[(l-2p)7-(l-4p)7±(2p)7],a2±-i(l- 
2p)2(l-4p)2[(l-2p)5±(l-4p)5],a3± = i(l-2p)4(l- 
4p)4[(l-2p)3±(l-4p)3],a4± = i(l-2p)6(l-4p)6[(l- 
2p) ± (1 — 4p)]. The coherent information can be quickly 
calculated. The typical time to obtain the optimal ?t.2 = 
133 is an hour (2004 CPU), while the former best result 
needs a week (2008 CPU)^. Further quick calculation 
is also possible by expanding the logarithmic in Eg. ipsj) . 
the time required for obtaining the optimal 712 for "7 in 
7^2" code is 2 minutes. For generic "77,1 in 712" code we 
have 



tCAT 



1 



nin2 



-{vq - Wl +7;2 - 1), 



with 
Vq = 
Vl = 

V2 = 



7i2(ao+ log2 ao+ -t- ai+ log2 ai+), 
772!(l-ao+-ai+)"^"^°"^^ 



E 



Mo!mi!("-2 -Mo -^i) 



c(mo,Mi), 



1 °° 1 

ln2 -^ 21(21- 1 



[(772-1)721 



21 



where c(mo,Mi) = 5(ao>i+ + ai+ao+) log2[5(ao+«i+ + 
)], 77^0 = 1, wi = 1, 7«, = 2Ci^{2 < j < 



•'i+"o+ 



^77i/2 



[(77-2 — 1)/2J , 7^771/2 = C'771 (for even 7^1). New results 
of optimal 712 are list in Table 1, with Pmax being the crit- 
ical value of channel noise such that / (pmax) = 0,and 
for p < Pmax we have positive coherent information. The 
"13 in 77-2" code has an optimal 712 > 1020. Too large 772 
makes the storage of Cii^^ overflow. 

The density matrix of the joint output can also be 
block diagonalizable for a non-additive code input which 
is diagonal in graph state basis. For example, the eigen- 
values of the output and joint output of ((5,6,2)) code 
with respect to depolarizing channel can be obtained ana- 
lytically when the input is diagonal and with equal prob- 
ability in graph state basis. The density matrix of the 
joint output can be block diagonalized as 32 blocks, each 
block is a 6 X 6 matrix and can be diagonalized eventually. 

In summary, we have block diagonalized the output 
density matrix of a code and the joint output density of 
the code and the ancilla system with respect to Pauli 
channel when the input quantum code is diagonal in 
graph state basis. For a {{n,L,d)) code which is diag- 
onal in graph state basis and encoding L states into 77, 
qubits with distance d, the joint output density matrix 
is reduced to 2" blocks, each is a L x L matrix. For a 
stabilizer code [[77, 1, d]] input which is diagonal in graph 



state basis and with an equal prior probability for all 2' 
encoded states, each block of joint output density matrix 
of a depolarizing channel can be further diagonalized with 
the Hadamard matrix. The eigenvalues are obtained in 
closed form, "ni in 712" concatenated repetition codes 
are used to illustrate the details. 
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